Phase-Ordering Dynamics with an Order-Parameter-Dependent Mobility: The 

Large-n Limit 
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The effect of an order-parameter dependent mobility (or kinetic coefficient), given by \(<j>) oc 
(1 — 4> 2 ) a , on the phase-ordering dynamics of a system described by an n-component vector order- 
parameter is addressed at zero temperature in the large-n limit. In this limit the system is exactly 
soluble for both conserved and non-conserved order parameter; in the non-conserved case the scaling 
form for the correlation function and it's Fourier transform, the structure factor, is established, with 
the characteristic length growing as L ~ ^/2(i+«)_ in tne conS e r ved case, the structure factor is 
evaluated and found to exhibit a multi-scaling behaviour, with two growing length scales differing 
by a logarithmic factor: Li ~ t 1/2( - 2+a) and L 2 ~ (t/ In t) 1/2( - 2+cl) . 
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I. INTRODUCTION 



In this paper we examine the effect of an order- 
parameter-dependent mobility, or kinetic coefficient, on 
the phase-ordering dynamics of a system described by an 
n-component vector order parameter. Both conserved 
and non-conserved order parameters are considered. For 
the case of a constant (i.e. order-parameter independent) 
mobility/kinetic coefficient, both these systems become 
analytically soluble in the large-n limit it is in 

this limit that we now consider the effect of an order- 
parameter-dependent mobility given by \(<fi) oc (1 — 4> 2 ) a , 
for models where the equilibrium order parameter satis- 
fies tf> = 1. Thus the mobility vanishes in equilibrium, 
leading to a reduction in the growth rate of the charac- 
teristic length scale, L(t), of the bulk phases. 

The effect of an order-parameter-dependent diffusion 
coefficient on a system with a scalar order-parameter has 
been studied by several authors || |(| since it has been 
proposed that for a scalar order-parameter a mobility of 
the form X((j)) = (1 — <f> 2 ) is required to accurately model 
the dynamics of deep quenches and the effect of an 
external field ||. Lacasta et. al. || studied this system 
numerically using a mobility given by A(<^>) = (1 — acf) 2 ). 
They found that for a = 1 the characteristic length grows 
as i 1 / 4 (in contrast to the conventional t 1 ^ 3 growth for 
a = 0), and for all a ^ 1 there is a crossover between 
L ~ t 1 / 4 and L ~ t 1 / 3 . Similar behavior was observed 
by Puri et al ||. This system has been solved exactly in 
the Lifshitz-Slyosov limit || for a more general mobility 
given by A(</>) = (1 — (f> 2 ) a ; in this system the system 
coarsens with growth exponent 1/(3 + a), despite the ab- 
sence of surface diffusion as a coarsening mechanism at 
late times (due to the geometry of the system), and the 
vanishing of the mobility in the bulk phases. 

Although a system described by a vector order param- 
eter will have a completely different morphology form 
the scalar case (e.g. there are no localized defects for 



n > d + 1), it is natural to try to generalise this order- 
parameter dependent mobility to the vector case 0. In 
this paper therefore we examine (in sections || and III) 
the coarsening dynamics of an n-component vector order- 
parameter for a general class of mobilities/kinetic coef- 
ficients given by \(<f>) = (1 — </> 2 ) Q , where a e for 
both the non-conserved and conserved cases. While these 
0(n) models are not exactly soluble for general n, exact 
solutions can be obtained in the limit n — > oo. 

In section || we consider a non-conserved system with 
a vector order parameter. The scaling hypothesis is es- 
tablished, and the exact forms of the two-time correlation 
function and the structure factor are calculated. We find 
that the characteristic length grows as L ~ i 1 / 2 ( 1 + Q ). 
Due to the absence of defects there is no Porod's law: 
the structure factor is Gaussian for all a. 



In the conserved case (section III), the structure fac- 
tor is found to depend on two characteristic lengths, 
Ll _ t i/2{2+ a ) and La _ (i/lnt) 1 /2(2+a) ) through the 



form S(k, t) - L 



d<t>(kL 2 ) 



This type of behaviour is termed 



'multiscaling', and the results for L\ and Li are general- 
izations of similar expressions obtained by Coniglio and 
Zannetti for the case of a constant mobility. Indeed, 
as expected, all the results of this paper reduce to the 
established constant A results when a is set to zero. 

We conclude with a summary and discussion of the 
results. 



II. THE NON-CONSERVED O(N) MODEL 

The dynamics of a non-conserved vector order pa- 
rameter are described by the phenomenological time- 
dependent Ginzburg-Landau equation [n0| , 



SFl 



A(<?) ( V'fr 



dV{<}> 2 ) 



(1) 



1 



where V(4> 2 ) is the potential energy term in the 
Ginzburg-Landau free energy functional, and is invari- 
ant under global rotations of (j>. In the following calcula- 
tion, the conventional choice is made for the form of the 
potential: 



(2) 



and the order-parameter-dependent kinetic coefficient is 
given by \{$) = (1 - 4> 2 ) a . 

In the limit n — > oo equation (Q) may be simplified by 
making the following substitution, 



-Hs.it 4) 



= n <cj>l>=(P) 



(3) 



where < ... > represents an ensemble average. Defining 
a(t) by the equation a(t) — (1 — n <4> 2 >), equation (|l|) 
then reduces to 



dt 



= a a {t)[V 2 + a(t) Ui. 



(4) 



If we now take the Fourier transform, this equation can 
easily be solved to give 



4 1) (t) = 4 l) (0)exp(-fc 2 6(t) + c(t)), 



(5) 



where b(t) = J*dt'a a (t') and c(t) = /„ dt' a 1+a {t'). On 
substituting equation (||) back into the definition of a(t) 
we find 



a(t) = 1 - Acxp[2c(£)] V exp[— 2k 2 b(t)], 



k 



(6) 



where we have used the conventional choice for the initial 
conditions. 



<4U%>-- 



A 



Sij k' 



(7) 



Using the fact that £ k exp(-2fc 2 6(i)) = {8irb(t))- d ' 2 in 
equation ([}]) we obtain 

o(t) = 1 - A[87r6(<)]- d / 2 exp[2c(i)]. (8) 

Since we are mainly interested in late times, we now solve 
this equation self-consistently to obtain the large-t result 
for b(t) and c(t). In order to make progress we make the 
assumption that at late times a(t) <C 1, and hence the 
term on the left-hand side of equation (ph may be ne- 
glected. The validity of this assumption will be proved a 
posteriori. Thus we wish to solve 



A[8vr&(i)]~ ex p[2c(t)] = 1. 



(9) 



Differentiating this expression with respect to time gives 
the following relation, 



%t) = 



d b{t) 
Jbjtj' 



(10) 



Substituting the derivatives of b(t) and c(i), which are 
given by: 



b(t)=a a (t), 
c(t)=a}+ a {t), 

into equation (|l0|), we find that 

, / s d 



(ii) 

(12) 



(13) 



If we now differentiate again, we obtain a simple differ- 
ential equation for a(t), and from this we find that the 
large-i behaviour of a(t) is given by 



a{t) 



4(l + a)A 1+ 



(14) 



Hence it can clearly be seen that the assumption that 
a(t) <C 1 at late times is justified. 

Using this result together with equations (^|) and (|l^) , 
we find that 



/ \ d 
c(t) 



4(1 + a) \t J' 



In I - 



where 



cr = (1 + a) !+° 



to 



d\~ 



a + l\ d 



A' 



1 \ l+a 



(15) 
(16) 



(17) 
(18) 



We are now in a position to evaluate the expression 
for the Fourier transform of the ord er p arameter at large 
t. Substituting equations ( |l5|) and (|16| ) into equation (g) 
we find that 



4 i) (t) = 4 i) (0)(^) 41+a exp(-<7^ 



(19) 



Using this result, we can evaluate the two-time structure 
factor and the correlation function. These are given by: 



S(k,t 1 ,t 2 



C(r,t u t 2 ) 



(8Tra) d ^(t 1 t 2 ) d /< 1+ ^ 
xexp(- ( Tfc 2 (^ /(1+Q) +4 /(1+a) )), (20) 

/ 4(t 1 ^ 2 ) 1 /( 1 +") 
\{t\ /{l+a) +t\ /{l+a) ) 2 ) 

: I , (21) 
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which, in the equal time case, reduce to the following 
expressions: 

S(k, t) = (Snaf^ 2 ^ exp (-2ak 2 t^ 1+a ^ , (22) 



C(r,t) = exp 



8at i/(i+a) 



(23) 



These results exhibit the expected scaling forms, with 
the characteristic length scale growing as L ~ ^ 1 / 2 ( 1 + Q ). 
The structure factor has a Gaussian form, without the 
power-law tail predicted by Porod's law. This is a direct 
consequence of the absence of defects in the system. 

If we now look at the two-time correlation function in 
the limit t\ 3> £2 1 we find that 



C{v,t u h) = 



"I 



!/(!+«)■ 



d/4 



exp 



4crq 



1/(1+") 



(24) 



Comparing this with the scaling form |l(J C(r, ti,fa) = 
(L2/ Li) x h(r / Li), we obtain the result, A = d/2, inde- 
pendent of a. 

It is also interesting to compare the response function, 
G(k,i) = (d4£\t)/d4$(Q)}, with the structure factor 

5(k, t, 0), i.e. with the correlation of <fi^ (t) with its i = 
value. Using equation (|l9| ) we find that: 



• l< " A(-J exp 



(-crfc 2 ^) , (25) 
(26) 



/ + \ 4(1 + Q ) / 

G(k,i) = - exp ^- CT fc 2 i— 



which verifies the relation S(k, t, 0) = AG(k, i). Note 
that this is an exact result valid beyond the large- n limit; 
this may be proved by integration by parts on the Gaus- 
sian distribution for {0k (0)} Jll| . 



III. THE CONSERVED O(N) MODEL 

The dynamics of a system described by a conserved 
vector order parameter are modelled by the Cahn- 
Hilliard equation ||l0f , 



= V. f A(<?)V ( -v 2 <^- 




(27) 



where we make the same choice for the potential as be- 
fore, V{4> 2 ) = |(1 — 4> 2 ) 2 . Following the method of the 

previous calculation, (j) 2 is eliminated using equation (|^) , 
therefore equation d27n reduces to 



dt 



= -a a (t)[V A & + a{t)V 2 



(28) 



where a(t) is defined as before. Taking the Fourier trans- 
form and solving the resulting differential equation yields 



4$ (*) = (0) exp (-k%(t) + k 2 c{t)) , 



(29) 



where b(t) and c(t) are defined as for the non-conserved 
case. Substituting this back into the formula for a(t) and 
using the random initial conditions given by equation ([?]) 
gives 

a(t) = 1 - A ^ exp (-2k 4 b(t) + 2k 2 c(t)) . (30) 



To make further progress we again assume that at large 
t, a(t) <C 1. This is checked for self-consistency later in 
the calculation. The sum over k is converted to an inte- 
gral and, using the change of variables 



equation (|30|) becomes 
A 



c{t) 



(31) 



where 



2 d - 1 n d / 2 T{d/2) \ b(t) J 

dxx^exp {2(3(t){x 2 -x 4 )) , 

0(t) = c 2 (t)/b(t). 



(32) 



(33) 



We now make an additional assumption (also to be 
verified a posteriori) that (3(t) — > 00 as t — > 00; the in- 
tegral on the left-hand side of equation ( |32| ) can then be 
evaluated by the method of steepest descents. Therefore, 
equation (B2h finally simplifies to 



A/3(t) 



-1/2 



23d/2 7r (d-l)/2l\ ( V2) \ b (t) J 



( ?W\ ' exp[/3(i)/2] = 1. 



(34) 



We now solve this equation asymptotically, obtaining 
expressions for a(t), b(t) and f3(t) at late times. On tak- 
ing the logarithm of equation (|34|) we find that 



/?(i)~-ln&(*) + 



2-d 



ln[ln &(£)]. 



(35) 



Using the definition of f3(t) (equation (p3|)) in equation 
( p5| ) we obtain an equation for c(t), which when differen- 
tiated, gives (to leading order) 



e(t) 



dhxb(t)Y /2 
8b{t) 



b(t). 



(36) 
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If we now substitute for the derivatives of b(t) and c(t) 
from equations (|ll| ) and (12) respectively, we find that 



a a (t) = b{t) 



dlnb(t)\ a/2 



8b{t) 

which has the asymptotic solution 



(37) 



b(t) 



(2 + a)t 



2/(2+a) 



dint 
4(2 + a) 



a/(2+a) 



(38) 



If we now differentiate this expression once more, we 
obtain the asymptotic behaviour of a(t), 



«(*) 



d\nt 



2{2 + a)H 



1/(2+") 



1 



2\nt 



and clearly a(t) C 1 at late times, justifying one of our 
initial assumptions. 

On substituting equation ( p8| ) into equation (|35|), we 
obtain 



lnt 



ln(mi). 



(40) 



We see that as t — > oo, fl(t) — > oo, justifying the appli- 
cation of the method of steepest descents to the integral 
in equation (p2[). Thus both our initial assumptions are 
satisfied. 

We are now in a position to evaluate the expression 
for 0^ (t) . Completing the square in the exponent on 
the right-hand side of equation (^9|) gives, 



^ ) W = ^ ) (0)exp 



m 

4 



■ 4 2 {(3(t) 



1/2 



(41) 



Substituting for b(t) and (3(t), from equations ( |38| ) and 
( fi"o| ) respectively, gives 



~ 4>^(0) (lni)^T t3(3TS) 



<p(k/k„ 



where 



h — 



dint 



2(2 + a) 2 t 



(42) 



(43) 



is the position of the maximum in the structure factor, 



2\2 



and 4>(x) = 1 — (1 — x 2 ) 

The structure factor is therefore given by 



S(k,t) ~ A(lni) 2 < 2 +°) t 2 < 2 ^ 



T<t>(k/k m ) 



(44) 



^From this expression it is self-evident that the struc- 
ture factor does not have the conventional scaling form 
S(k, t) ~ L d g(kL). In this system there are two different 



length scales, L\ and L 2 , which which differ only by a 
logarithmic factor and are given by 



Ll „ t l/2(2+a) > 



k^ 1 - 



t 

hit 



l/2(2+a) 



(45) 
(46) 



The structure factor is therefore of the form 5(k, i) 



with an additional logarithmic correction fac- 

2 + a-d 

tor, (lnt) 2 ( 2 +°) ; the exponent depends continuously on a 
scaling variable. This type of behaviour is called 'multi- 
scaling', and was first noted by Coniglio and Zannctti for 
the case a = Rl. Note that the a-dependence enters 
through the length scales Li and L 2l while the function 
(39) cf)(x) is independent of a. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper we have considered the effect of an order- 
parameter-dependent mobility/kinetic coefficient, given 
by A(</>) = (1 — <p 2 ) a , on a system described by an n- 
component vector order parameter. Exact results have 
been obtained in the large- n limit, a limit which de- 
spite its limited applicability to physical systems has been 
widely studied as one of the few exactly soluble models 
of phase-ordering kinetics [^]-^|Jl^-|l^1 . All the results ob- 
tained reduce to the expected constant A results when a 
is set to zero. 

In the non-conserved system, the correlation function 
and its Fourier transform, the structure factor, were 
explicitly calculated and found to be of the expected 
scaling form, with the characteristic length growing as 
L ~ £i/ 2 (i+"). The order-parameter-dependent kinetic 
coefficient slows down the rate of domain coarsening; the 
result reduces to the familiar i 1 / 2 growth for the case 
a = The result A = d/2, independent of 

a, was established from the two-time correlation func- 
tion C(r, ttjti) in the regime t\ 3> £2, and the relation 
S(k, t,0) — AG(k, t), relating the correlation with, and 
the response to, the initial condition was verified. The 
equal-time correlation functions and structure factor are 
Gaussian. 

The system with a conserved order parameter was 
found to exhibit a more unusual behaviour. In this sys- 
tem, the structure factor does not have the conventional 
scaling form and is dependent on two scaling lengths, 

t l/2(2+a) and k -l _ (t/lnt)l/2(2+«) ) where km is thc 

position of the maximum in the structure factor. This 
type of behaviour was first discovered in a phase-ordering 
system by Coniglio and Zannetti 0, for the a = 
case. For a = this behaviour is a consequence of 
the non-commutativity of the large-n and large-t limits, 
as demonstrated within a soluble approximate model by 
Bray and Humayun [Q. They demonstrated that for fi- 
nite n, in thc limit t — ■> 00, conventional scaling is found 
whereas if the n — > 00 limit is taken first (at finite t), the 
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Coniglio and Zannetti result |jj is recovered. At large, 
but finite n, multiscaling behaviour is found at interme- 
diate times, with a crossover to simple scaling behaviour 
occurring at late times jl4|,|l7],[l8) . We anticipate that a 
similar crossover to simple scaling at late times will occur 
for any a for large but finite n, leaving a single growing 
length scale L ~ i 1 / 2 ( 2 +a) j j-,^ an explicit demonstration 
of this goes beyond the scope of the present work. 

Note that all the results presented above have been de- 
rived in the absence of thermal noise, so these results are 
strictly valid only for quenches to T = 0. However, since 
we do not expect temperature to be a relevant variable 
jl0|,[l2]], qualitatively similar results should be obtained 
for quenches to T > (but T < T c ), at least for noncon- 
served dynamics (with n finite or infinite) or conserved 
dynamics with finite n |n| . 
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